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Abstract: We show how the infinite color-Coulomb energy of color-charged states is re- 
lated to enhanced density of near-zero modes of the Faddeev-Popov operator, and calculate 
this density numerically for both pure Yang-Mills and gauge-Higgs systems at zero temper- 
ature, and for pure gauge theory in the deconfined phase. We find that the enhancement 
of the eigenvalue density is tied to the presence of percolating center vortex configurations, 
and that this property disappears when center vortices are either removed from the lattice 
configurations, or cease to percolate. We further demonstrate that thin center vortices have 
a special geometrical status in gauge-field configuration space: Thin vortices are located at 
conical or wedge singularities on the Gribov horizon. We show that the Gribov region is 
itself a convex manifold in lattice configuration space. The Coulomb gauge condition also 
has a special status; it is shown to be an attractive fixed point of a more general gauge 
condition, interpolating between the Coulomb and Landau gauges. 
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1. Introduction 



Of the many different ideas that have been advanced to explain quark confinement, more 
than one may be right, or at least partially right, and these should be related in some 
way. As in the old story of six blind men describing an elephant, those of us concerned 
with the QCD confinement mechanism might benefit from unifying some of our separate 
impressions, to form a better image of the entire beast. 

In this article we investigate the relationship between center vortices (for a review, 
cf. ref. [jij) and the Gribov horizon in Coulomb gauge, whose relevance to confinement 
has been advocated by Gribov and Zwanziger Q. We begin with the simple fact that in 
a confining theory, the energy of an isolated color charge is (infrared) infinite, and this 
energy is a lower bound on the color-Coulomb energy || (defined below). This fact implies 
an enhancement of the density of near-zero eigenvalues A n 

M<P (n) = \ n (j) {n) (1.1) 

of the Faddeev-Popov (F-P) operator in Coulomb gauge 

M = -V -V(A), (1.2) 

where T> (A) is the covariant derivative. The F-P eigenvalue density is an observable which 
we are able to calculate numerically (section ^), via lattice Monte Carlo and sparse-matrix 
techniques. Applying standard methods j| we are able to separate any thermalized lat- 
tice configuration into vortex-only and vortex-removed components, and we find that the 
enhancement of the F-P eigenvalue density can be entirely attributed to the vortex compo- 
nent of the gauge fields. Vortices are associated with an enhancement in the density of F-P 
eigenvalues A near A = 0; this enhancement is key to the divergence of the color-Coulomb 
energy in Coulomb gauge. It is absent in the Higgs phase of a gauge-Higgs system (section 
|3|), where remnant gauge symmetry is broken [0], and vortices cease to percolate. In partic- 
ular, we compare the F-P eigenvalue density found in the Higgs phase to the corresponding 
density for configurations in the confined phase, with vortices removed. These densities are 
identical in form. The density of F-P eigenvalues in the high-temperature deconfined phase 
of pure gauge theory is examined in section |j. We find that the linearly rising, unscreened 
color-Coulomb potential, which is present in the deconfined phase ||, is associated with 
an enhanced density of F-P eigenvalues, and that this enhancement is again attributable 
to the vortex component of the gauge field. Divergence of the color-Coulomb energy is a 
necessary but not sufficient condition for color confinement, and this phase provides an in- 
teresting example where the infinite color-Coulomb potential gets screened. Although the 
horizon scenario was invented to describe confinement, it nicely accounts for the divergence 
of the color-Coulomb energy in the deconfined phase, as explained in the Conclusion. It is 
also shown (section ||) how an array of center vortices, set up "by hand" to simulate some 
aspects of a percolating vortex configuration, leads to an accumulation of F-P eigenvalues 
near A = 0. 

In section ^ we demonstrate that center configurations (equivalently: thin center vor- 
tices) have some remarkable geometrical properties in lattice configuration space. First, as 
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already shown in ref. ||, center configurations lie on the Gribov horizon. It is known that 
the Gribov horizon is a convex manifold for continuum gauge fields || — a result which 
we extend here (section ||) to lattice gauge fields — and one might therefore suspect that 
the Gribov horizon is also smooth and differentiable. In fact, thin vortex configurations 
turn out to be distinguished points on the Gribov horizon, where the manifold acquires 
conical or "wedge" singularities. Finally, in section ||, we point out that the Coulomb gauge 
condition also has a special status, in that Coulomb gauge is an attractive fixed point of 
a more general interpolating gauge condition, which has the Coulomb and Landau gauge 
conditions as special cases. 



2. The F-P eigenvalue density 

The energy of an isolated color charge is (infrared) infinite in the confined phase, even on 
the lattice where ultraviolet divergences are regulated. We will consider charged states in 
Coulomb gauge, which, for a single static point charge, has the simple form 

*%[A;x]=i> a {x)MA\, (2.1) 

where a is the color index for a point charge in color group representation r, and is the 
Coulomb gauge ground state. The gauge-field excitation energy £ r of this state, above the 
ground state energy, is due entirely to the non-local part of the hamiltonian 

£r ~ (*Z[A ]X ]\*Z[A;x]) ~ ^Hcout^o). (2.2) 
We recall that in Coulomb gauge, the hamiltonian is a sum H = H g i ue + H cou i, where 
H g i U e = ^Jd 3 x (J-^E^ a J ■ fi**J-\ + B a ■ B a ), 
ticoui = \ f d 3 xd 3 yJ-hp a (x)JK ab (x,y;A)p b (y)J-^ 



and 



K ab (x,y;A) = [M _1 (-V 2 )Ar 



r -ll ab 

J xy ' 

„a a _ fabc Ab zritr.c 

P ~ Pq~ 9 J A k h k , 

J = det[-V • T>{A)]. (2.3) 



Then the excitation energy of the charged state, in SU(iV) gauge theory, is 

£ r =g 2 ±-Tr[L a L b ](K ab (x,x;A)) = ^^^-(^(x.x; A)), (2.4) 

where the {L a } are color group generators, d r the dimension, and C r the quadratic Casimir 
of representation r of the color charge. Therefore, the excitation energy £ r is the energy of 
the longitudinal color electric field due to the static source, which we can identify as the 
color Coulomb self-energy. This energy is ultraviolet divergent in the continuum, but of 
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course that divergence can be regulated with a lattice cut-off. The more interesting point 
is that in a confining theory, £ r must still be divergent at infinite volume, even after lattice 
regularization, due to infrared effects. 

The excitation energy ( |2.2| ) represents the (infrared-infinite) energy of unscreened color 
charge in the state fl2,l| ), which in general is a highly excited state that is not an eigenstate 
of the hamiltonian. States of this kind are useful for extracting the self-energy of an 
isolated charge due to its associated color-Coulomb field. On the other hand, the minimal 
free energy £ s of a state containing a static external charge, at inverse temperature T, 
is obtained from the value of the Polyakov loop P ~ exp(— £ S T). This minimal energy 
may be infrared finite in an unconfined phase even if £ r is not, providing the external 
charge can be screened by dynamical matter fields, or by high-temperature effects. The 
infrared divergence of £ r must be understood as a necessary, but not sufficient, condition 
for confinement. 

We note in passing that the charged state ( |2.1[ ) in Coulomb gauge corresponds, in QED 
in temporal gauge, to the well-known form 



V QED [A;x] =exp 



[ d 3 z A(z) ■ V— 
J 4tt\x - 



j>(x)*3 ED [A] (2.5) 



The investigation of this type of "stringless" state with external charges in non-abelian 
theories, using perturbative methods, was undertaken some time ago by Lavelle and Mc- 
Mullan in ref. 0. The exponential prefactor in eq. (|2.5[) can be identified as the gauge 
transformation taking an arbitrary configuration A k {x) into Coulomb gauge. This feature 
generalizes to non-abelian theories, and "stringless" states with static charges in temporal 
gauge can be formally expressed in terms of the gauge transformation to Coulomb gauge, 
as shown in ref. |J. 

We now proceed to the lattice formulation, with an SU(2) gauge group. The link 
variables can be expressed as 

U„(x) = b^x) + iff ■ a^x) , b^xf + Y,al{xf = 1 (2.6) 

c 

and (when the lattice version of the Coulomb gauge condition V • A = is satisfied) the 
lattice Faddeev-Popov operator is 

3 

M± = 5 ab Y, {Sxy [bk(x) + b k {x - k)] - S^hix) - 8 y x _ k b k (y)} 



k=l 

3 

„abc 



£{^4(*)-^_*4(s/)}, (2-7) 



fc=l 

where indices x,y denote lattice sites at fixed time. Denote the Green's function corre- 
sponding to the F-P operator as 1 

a(n) ib(n)* 



1 This expression assumes M is invertible, a point which will be discussed in subsection 2.1 



below. 
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where (ftx^ , A ra are the n-th normalized eigenstate and eigenvalue of the lattice F-P oper- 
ator M£y. Defining the representation-independent factor £ = £ r /(g 2 C r ) in the Coulomb 
self-energy, we find 

■(K aa (x,x;U)) 



N 2 

= o y~^(^?/i(~'^ 7 )yiy2^y2x) = 7yy~ ^ ^ (G xyi (— V }yxy-fiy\x) 
yiy-2 3 x 2/12/2 

1 I Mm) Mm)* Mn)Mn)*\ 

_ 1 Y^ Y^ Y^ Y^ / ijE gyi <■ v ) y2 \ 

- 3y3 2^2^2^\ Am An / 

6 x yiy 2 rn n \ m 71 / 

_ 1 Y^ Y^ / fyi l~ v )mm9y2 \ / 9 q\ 

~W^^\ aI /' ( '- 9j 

2/12/2 n \ " / 

where V3 = L 3 is the lattice 3-volume. Also defining 

F n = Y. < n) *(- v2 ).i. 2 < n) = fi n> ■ (-V 2 )<? (n) (2.10) 

3/13/2 

we have 



n 



The Faddeev-Popov operator, on the lattice, is a 3V3 x 3V3 sparse matrix; the number 
of linearly-independent eigenstates is therefore 3V3. Let iV(A,A + AA) be the number of 
eigenvalues in the range [A, A + AA]. We define, on a large lattice, the normalized density 
of eigenvalues 

, , N(X, A + AA) 

" (A) s 3V3AA • < 2 - 12) 

Then as the lattice volume tends to infinity, 

£ = J o ^ p(A)*m (2.13) 

where it is understood that the integrand is averaged over the ensemble of configurations. 
From this we derive a condition for the confinement phase: The excitation energy £ r = 
g 2 C r £ of a static, unscreened color-charge is divergent if, at infinite volume, 

Dm " W f (A > > 0. (2.14) 
A^o A 

In perturbation theory, at zero-th order in the gauge-coupling, the Faddeev-Popov operator 
is simply a laplacian, whose eigenstates are plane waves. Then A = k 2 , where k is the 
momentum, and from the volume element of momentum space it is easy to see that to 
zeroth order, in the limit of infinite lattice volume, 

p(X) = ^, F(A) = A, (2.15) 
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which obviously does not satisfy the confinement condition ( 2.14 ). At zeroth-order we have 
£ = Amax/(2vr 2 ). 

We now recall briefly some aspects of the Gribov horizon scenario Q . The lattice ver- 
sion of the Coulomb gauge-fixing condition V • A = is satisfied by any lattice configuration 
U such that 

3 

R = Y,Y, TY ^{x)U k {x)g-\x + k)] (2.16) 

x k=l 

is stationary at the trivial gauge transformation g = I. The Faddeev-Popov operator is 
obtained from the second derivative of R with respect to gauge transformations, so if R is a 
local maximum, then all the eigenvalues of the Faddeev-Popov operator are positive. The 
subspace of configuration space satisfying this condition is known as the Gribov region, 
and it is bounded by the Gribov horizon, where M develops a zero eigenvalue. In principle, 
the functional integral in minimal Coulomb gauge should be restricted to a subspace of the 
Gribov region in which the gauge fields are global maxima of R; this subspace is known 
as the "Fundamental Modular Region." Part of the boundary of the fundamental modular 
region lies on the Gribov horizon. 

The dimension of lattice configuration space is of course very large, on the order of the 
number of lattice sites, and it has been proven that, in contrast to an abelian theory, the 
Gribov region is bounded and convex (cf. Q and section || below). A suggestive analogy is 
that of a sphere in iV-dimensions, which has a volume element proportional to r ^dr, so 
that most of the volume is concentrated close to the surface of the sphere. By this analogy, 
it is reasonable to suppose that the volume of the Gribov region is concentrated near the 
Gribov horizon. If that is so, then typical configurations generated by lattice Monte Carlo, 
fixed to Coulomb gauge by standard over-relaxation techniques, will also lie very close to 
the Gribov horizon, and the F-P operator for such configurations will acquire near-zero 
eigenvalues. This is not enough by itself to ensure confinement; even the laplacian oper- 
ator will have a spectrum of near-zero modes at large lattice volumes. The conjecture is 
that typical configurations near the Gribov horizon may also have enhanced values for p(A) 
and F(X) (compared to the perturbative expressions) at A — > 0, such that the confinement 
condition ( [2.14D is satisfied. Our task is to check, numerically, whether this enhancement 
exists or not. 



2.1 Observables 

We apply the ARPACK routine || , which employs a variant of the Arnoldi procedure for 
sparse matrices, to evaluate the lowest 200 eigenvectors and corresponding eigenvalues of 
the F-P matrix Mf y in eq. (|2.7] ), for configurations generated by lattice Monte Carlo. The 
first three eigenvalues are zero for SU(2) gauge theory, regardless of the lattice configu- 
ration, due to the fact that the eigenvector equation ( |1 . 1[) is trivially satisfied by three 
linearly independent, spatially constant eigenvectors 

^) = ^=5 an , n = l,2,3 (2.17) 
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with zero eigenvalue. 

The existence of these trivial zero modes is related to the fact that physical states 
with non-zero total color charge cannot exist in a finite volume with periodic boundary 
conditions. This is true even for an abelian theory, and the reason is simple: the Gauss 
law cannot be satisfied for total non-zero charge in a finite periodic lattice. In such cases, 
the electric flux lines diverging from point sources have nowhere to end. This means that 
the F-P operator (or the laplacian, in an abelian theory) is non-invertible on a periodic 
lattice. It is precisely the existence of the trivial zero modes which makes the F-P operator 
non-invertible; there is no such difficulty in an infinite volume. In order to extrapolate 
our results on finite volumes to infinite volume, which allows non-zero total charge, there 
are two possibilities. First, we could get rid of zero modes by imposing non-periodic (e.g. 
Dirichlet) boundary conditions on the finite lattice. Second, we could perform our Monte 
Carlo simulations on finite periodic lattices as usual, but drop the trivial zero modes before 
extrapolating our results to infinite volume. In this article we choose the second approach, 
and exclude the trivial zero modes from all sums over eigenstates. 

The average eigenvalue density p(X) is obtained from the remaining 197 eigenvalues in 
each thermalized configuration (there are L such configurations in a given L 4 lattice, one at 
each time-slice). The range of eigenvalues is divided into a large number of subintervals, and 
eigenvalues are binned in each subinterval to determine the expected number N(X, A + AA) 
of eigenvalues per configuration falling into each bin. Substituting this value into the 



definition ( 2.12 ) of the normalized density of states, we obtain an approximation to p(X) 
for A values in the middle of each subinterval. We also compute the expectation value of 
the n-th eigenvalue and corresponding quantity -F(A n ) 

<A n > and (F n ) = • (-V 2 )<? (n) ) (2.18) 

for the n = 4 — 200 non-trivial eigenvectors. Our plots of F(X) vs. A, shown below, are 
obtained by plotting (F(X n )) vs. (A n ). Finally, we calculate the average contribution of 
low-lying eigenstates with A ra < 0.15 to the energy £ of unscreened color charge: 

For our purposes, the precise value of the upper limit in the sum is not too important. 
We have chosen the upper limit A = 0.15 in order that the 200 lowest eigenvalues, on each 
lattice volume we have studied, include the range < A < 0.15. 

2.2 Center projection and vortex removal 



Not every configuration on the Gribov horizon satisfies the confinement condition ( 2.14| ). 



For example, any purely abelian configuration in a non-abelian theory lies on the Gribov 
horizon (after gauge transformation to Coulomb gauge ||) and therefore has a (non-trivial) 
zero eigenvalue, but not all such configurations will disorder Wilson loops, or lead to an 
F-P eigenvalue spectrum satisfying eq. ( |2.14j ). The center vortex theory of confinement 
holds that a particular class of field configurations dominates the vacuum state at large 
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scales, and is responsible for the linear rise of the static quark potential. If so, then these 
same configurations should be responsible for the pileup of F-P eigenvalues near A = 0, 
resulting in the infinite energy of an isolated color charge. This is the connection which we 
think must exist between the center vortex and Gribov horizon confinement scenarios. 

To investigate this connection, we apply standard methods to factor a lattice configu- 
ration into its vortex and non- vortex content. This is done by first fixing to direct maximal 
center gauge, which is the Landau gauge condition in the adjoint representation 

R = ^2 T* - ^^)] 2 is maximum (2.20) 

X fl 

using an over-relaxation technique. The lattice configuration is factored into 

U^x) = Z^xfi^x) (2.21) 

where 

Z^x) = sign{Tr[[/ M (x)]} (2.22) 

is the center-projected configuration, and f7„ is the "vortex-removed" configuration. The 
center-projected (thin vortex) configuration Z^(x) carries the fluctuations which give rise 
to an area law for Wilson loops. The asymptotic string tension of the vortex-removed 
configuration U^x) vanishes, as does its chiral condensate and topological charge. The 
numerical evidence supporting these statements is reviewed in ref . |I| . 

In our procedure, each thermalized lattice configuration is transformed to direct max- 
imal center gauge and factored as above into a center-projected configuration Z^(x), and a 
vortex-removed configuration U^(x). These are then transformed separately into Coulomb 
gauge. Of course, any center-projected configuration Z /1 (x), with links = ±7, trivially 
fulfills the Coulomb gauge condition 

^Tr[o- a (U k (x) + ul(x-k))j =0, (2.23) 

k 

but in general such configurations are far from the minimal Coulomb gauge, and are nor- 
mally not even in the Gribov region. So in practice we perform a random gauge trans- 
formation on Z^(x), and then fix to a gauge copy in the Gribov region by the usual 
over-relaxation method. We will refer to such copies as "vortex- only" configurations. Ap- 
plying the Arnoldi algorithm to calculate the F-P eigenvectors and eigenvalues, we compute 
observables e, p(A), (F n ), (X n ) for both the vortex-only and vortex-removed configurations. 

Any purely center configuration lies on the Gribov horizon. By "center configuration" 
we mean a lattice all of whose link variables can be transformed, by some gauge transfor- 
mation, to center elements of the gauge group, and the "vortex-only" configurations are 
center configurations in this sense. It was shown in ref. ||] that for the SU(2) gauge group, 
such configurations have in general three non-trivial F-P zero modes, in addition to the 
three trivial, spatially constant zero modes ( 2. 17] ), 2 In computing e for the vortex-only 



configurations, we therefore exclude the first six eigenvalues. 



2 In some special cases of high symmetry there may be fewer non-trivial zero modes. 
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2.3 Results 



Most of our simulations have been carried out at (3 = 2.1. This is not such a weak 
coupling, but it allows us to use modestly sized lattices whose volumes are nonetheless 
quite large compared to the confinement scale, and to study the volume dependence. The 
F-P observables are calculated in the full configurations, the thin-vortex configurations, 
and the vortex-removed configurations, each of which has been transformed to Coulomb 
gauge. 

Figures || and § show our results for p{\) and F(X) for the full configurations, on a 
variety of lattice volumes ranging from 8 4 to 20 4 (to reduce symbol overlap near A = 0, 
we do not display the entire set of data points in .F(A)). The apparent sharp "bend" in 
p{\) near A = becomes increasingly sharp, and happens ever nearer A = 0, as the lattice 
volume increase. The impression these graphs convey is that in the limit of infinite volume, 
both p(A) and F(X) go to positive constants as A — > 0. However, for both p(A) and -F(A) 
we cannot exclude the possibility that the curves behave like A p , X q near A = 0, with p, q 
small powers. If we assume that the low-lying eigenvalue distribution scales with the total 
number of eigenvalues (3L 3 ) in the manner suggested by random matrix theory, then it is 
possible to deduce, from the probability distribution of the lowest non-zero eigenvalue, the 
power dependence of p{\) near A = 0. This analysis is carried out in Appendix |A], and 
gives us the estimates 

p{\) ~ A ' 25 , F(X) ~ A ' 38 full configurations (2.24) 

at small A and large volume, with perhaps a 20% error in the exponents. With this behavior 
the Coulomb confinement condition is satisfied, and the Coulomb self-energy is infrared 
divergent. 

In Fig. [H we plot e vs. lattice extension L, together with a best straight-line fit through 
the points at L = 10 — 20. The cut-off energy e rises with L, and this rise is consistent with 
linear, although the data is not really good enough to determine the precise L-dependence. 




Figure 1: The F-P eigenvalue density at f3 = 2.1, on 8 4 — 20 4 lattice volumes. 
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Figure 2: F(X), the diagonal matrix element of (—V 2 ) in F-P eigenstates, plotted vs. F-P 
eigenvalue. 

One might wonder how it is possible that a pair of quark-antiquark charges, in a global 
color singlet combination, can have a finite total Coulomb energy in an infinite volume, 
given that each charge has a divergent self-energy. This question was addressed in ref. Q, 
which computed Coulomb energies from timelike link correlators. The answer is that both 
the quark self-energies, and the energy due to instantaneous one-gluon exchange between 
separated sources, contain an infrared divergent constant. It can be shown that for global 
color singlets these constants precisely cancel, while in non-singlets the self-energy is not 
entirely cancelled, and the total energy is infrared divergent. 

Next we consider the F-P observables for the "vortex-only" configurations, consisting 
of thin vortex configurations (in Coulomb gauge) which were extracted from thermalized 
lattices as described above. Our data for p(X) and F(X) at the same range (8 — 20 4 ) of 



4 
3.5 
3 

w 2.5 
2 
1.5 
1 



full configurations 



10 15 20 25 
L 



Figure 3: e vs. lattice extension L, for full, unprojected lattice configurations. The straight 
line is a linear fit through points at L = 10 — 20. 
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Eigenvalue density, (5=2.1 




Figure 4: F-P eigenvalue density in vortex-only configurations. 

lattice volumes is displayed in Figs. ||| and ||. The same qualitative features seen for the full 
configurations, e.g. the sharp bend in the eigenvalue density near A = 0, becoming sharper 
with increasing volume, are present in the vortex-only data as well, and if anything are 
more pronounced. From the graphs, it would appear that 

p(0) 0.06 , F(0) » 1.0 , vortex-only. (2.25) 

This non-zero limit for /9(A), F(X) at A — > is supported by an analysis of the low- 
eigenvalue universal scaling behavior as a function of L, which is reported in Appendix [A|. 
Once again, the confinement criterion ( |2.14| ) is obviously satisfied. 

Figure || shows our data for e(L), again with a linear fit through the data points at 
L = 10 — 20, although the linear dependence is not really established. 




Figure 5: F(X), the diagonal matrix element of (—V 2 ) in F-P eigenstates, for vortex-only 
configurations. 



vortex only 
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Figure 6: e vs. lattice extension L, for the vortex-only lattice configurations. The straight 
line is a linear fit through points at L = 10 — 20. 



Finally, we consider the F-P observables of the vortex-removed configurations U trans- 
formed to Coulomb gauge. Our results are shown in Fig. ^ for p(X), and Fig. || for F(X). The 
behavior of these observables is strikingly different, in the vortex-removed configurations, 
from what is seen in the full and vortex-only configurations. A graph of the eigenvalue 
density, at each lattice volume, shows a set of distinct peaks, while the data for F(X) is 
organized into bands, with a slight gap between each band. Closer inspection shows that 
eigenvalue interval associated with each band in F(X) precisely matches the eigenvalue 
interval of one of the peaks in p(X). 

In order to understand these features, we con- 
sider the eigenvalue density of the F-P operator M°J = 
5 ab (—V 2 ) xy appropriate to an abelian theory (or a 
non-abelian theory at zero-th order in the coupling) . 
Although we can readily derive the result p(X) ~ \fX 



(2.15) at infinite volume, this result is slightly mis- 
leading at finite volume, where the eigenvalue den- 
sity is actually a sum of delta-functions 
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(1,1,1) 


12sin 2 (vr/L) 
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18 


(2,0,0) 


4sin 2 (27r/L) 



(2.26) 



Table 1: Eigenvalues of the zero- 
field lattice F-P operator -<5 a6 V 2 , 
and their degeneracies N^. 



In the sum, the index k labels the distinct eigenvalues of the lattice laplacian —V 2 , and 
is the degeneracy of Xk- Explicitly, to each distinct eigenvalue Xk on an L 4 lattice, there is 
a set of integers ni_3 such that 



sin 



i=l 



TTUi 

L 



(2.27) 



The first few values of Xk, and their degeneracies, are listed in Table [|. 
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vortex-removed configurations 



L=20 



0.05 0.1 0.15 0.2 0.25 

X 



Figure 7: F-P eigenvalue densities for vortex-removed configurations, on a 20 4 lattice 
volume. 

We have compared these degeneracies Nk, of the zeroth-order F-P eigenvalues, with the 
number of eigenvalues per lattice configuration found inside the fc-th "peak" of p(\), and k- 
th "band" of F(X). We find there is a precise match. This leads to a simple interpretation: 
the vortex-removed configuration Uu can be treated as a small perturbation of the zero-field 
limit Ufj, = I. This perturbation lifts the degeneracy of the spreading the degenerate 
eigenvalues into the bands of finite width in A. At least for small k, these bands do not 
overlap. Likewise, the perturbation broadens the infinitely narrow <5-function peaks in the 
density of eigenstates, eq. Q2.26D , into the peaks of finite width seen in Fig. ^. 

Because both the density of eigenvalues and the data for F(\) seem to be only a 
perturbation of the corresponding zero-field results, it appears to be most unlikely that the 
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Figure 8: .F(A) vs. A in the vortex-removed configurations. We display results for the 20 4 
volume alone (left), and a variety of lattice volumes (right). 
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no- vortex configurations lead to a divergent Coulomb self-energy. In fact, we find that the 
low-lying eigenvalue spectrum scales with lattice extension L as 

(An) ~ -Jj, (2.28) 

just as in the zero-field (or zero-th order) case. We have not plotted e(L) for the vortex 
removed case, because for the smallest lattice volume no eigenvalues actually lie in the 
given range < A < \ m i n , and even at L = 10, 12 only a few of the lowest eigenvalues are 
in this range. We do note, however, that the values of e on the large lattices are roughly 
two orders of magnitude smaller than the corresponding values for the full configurations. 

Vortex removal is in some sense a minimal disturbance of the original lattice configu- 
ration. The procedure only changes field strengths at the location of P-vortex plaquettes, 
and the fraction of P-plaquettes out of the total number of plaquettes on the lattice dives 
exponentially to zero with increasing f3. Nevertheless, this modest change of the lattice con- 
figuration is known to set the string tension and topological charge to zero, and to remove 
chiral symmetry breaking. At (3 = 2.1, we have found that vortex removal also drastically 
affects the density of low-lying eigenvalues, and one sees, in the multi-peak structure of 
/9(A), the remnants of the delta- function peaks of the free theory. We have used this com- 
paratively low value of (3 so that we may probe large physical distances, as compared to 
the distance scale set by the string tension, on rather modestly sized (up to 20 4 ) lattices. 
This allows us carry out the finite volume scaling analysis reported in Appendix A. 

However, using such a small (3 has a price, in terms of our confidence in the effects 
of vortex removal. At (3 = 2.1 roughly 17% of all plaquettes are P-vortex plaquettes (cf. 
Fig. 20 in ref. |l|]), and one may object that in this case vortex removal is not such a small 
disturbance of the lattice configuration. Perhaps the drastic effect of vortex removal on the 
eigenvalue density is simply an artifact of the substantial number of plaquettes modified. 

In order to address this concern, we have computed p(X) and -F(A) at (3 = 2.3 and 
j3 = 2.4, where the P-vortex densities have dropped to around 9% and 4%, respectively, 
of the total number of plaquettes. The results are shown in Fig. [9|, where we display the 
data for the unmodified, vortex-only, and vortex-removed configurations on the same plot. 
Only two lattice volumes are shown at each coupling. The effect of vortex removal is seen 
to be much the same at these higher (3 values as at (3 = 2.1. Again we see a multi-peak 
structure in p(X), and a band structure in -F(A) at the low-lying eigenvalues (although the 
gap between bands narrows as (3 increases). In each configuration, we have checked that the 
number of eigenvalues in each peak of p(X) matches the number of eigenvalues in each band 
of F(X), and that this number is again equal to the degeneracy of eigenvalues. Therefore, 
the interpretation proposed for the no-vortex data at /3 = 2.1 appears to apply equally 
well at the higher (3 values: these data are simply perturbations of the zero-field result, for 
which the eigenvalue density is a sum of delta-functions. The data for the vortex-only and 
unmodified configurations are also qualitatively similar to the results we have obtained at 
the lower (3 = 2.1 value, although the lattice volumes, in physical units, are considerably 
smaller than the corresponding lattice volumes at (3 = 2.1. 
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Figure 9: p(A), i^(A) at couplings f3 = 2.3, 2.4 for full, vortex-only, and vortex-removed 
configurations. 



We conclude that it is the vortex content of the thermalized configurations which is 
responsible for the enhancement of both the eigenvalue density and F{X) near A = 0, 
leading to an infrared-divergent Coulomb self-energy. 

3. Gauge-Higgs theory 

Next we consider a theory with the gauge field coupled to a scalar field of unit modulus in 
the fundamental representation of the gauge group. For the SU(2) gauge group, the lattice 
action can be written in the form || 

S = P^2h Tr P Uuju ^ + 7 Yl ^[^{^U^x^ix + V)] (3.1) 

plaq x,ix 

with (j) an SU(2) group-valued field. Strictly speaking, this theory is non-confining for all 
values of /3, 7. In particular, there is no thermodynamic phase transition (non-analyticity 
in the free energy), from the Higgs phase to a confinement phase; this is the content of a 
well-known theorem by Osterwalder and Seiler |jTo| l, and Fradkin and Shenker [11]. 
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However, the Osterwalder-Seiler — Fradkin-Shenker (OS-FS) theorem is not the last 
word on phase structure in the gauge-Higgs system. In fact, there is a symmetry-breaking 
transition in this theory. If one fixes to Coulomb gauge, there is still a remaining freedom 
to perform gauge transformations which depend only on the time coordinate. It is there- 
fore possible, in an infinite volume, that this symmetry is broken on any given time-slice 
(where the remnant symmetry is global). The order parameter for the symmetry-breaking 
transition is the modulus of the timelike links, averaged over any time slice 



V = pj>„(x,t). (3.2) 



If Q — > as L — > oo, then the remnant gauge symmetry is unbroken, the Coulomb potential 
(as opposed to the static potential) between quark-antiquark sources rises linearly, and the 
energy of an isolated color-charge state of form ( |2.1| ) is infrared divergent. Conversely, if 
Q > at infinite volume, then the remnant symmetry is broken, the Coulomb potential is 
asymptotically flat, and the energy of an isolated color-charge state is finite. 

These matters are explained in some detail in ref. Q, where we report on a sharp 
transition between a symmetric ("confinement-like") phase and a broken ("Higgs") phase 
of remnant gauge symmetry. There is no inconsistency with the OS-FS theorem, which 
assures analyticity of local, gauge-invariant order parameters. The order parameter Q, 
when expressed as a gauge-invariant observable, is highly non-local. The transition line 
between the symmetric and broken phases is also not the location of a thermodynamic 
transition, in the sense of locating non-analyticity in the free energy. Rather, it is most 
likely a Kertesz line, of the sort found in the Ising model at finite external magnetic field, 



which identifies a percolation transition [12|. In the gauge-Higgs case, the objects which 



percolate in the confinement-like phase, and cease to percolate in the Higgs phase, turn 
out to be center vortices |1| (see also [14, 15]). 



In the previous section we investigated the effect, on F-P observables, of removing 
center vortices from lattice configurations by hand. The gauge-Higgs system gives us the 
opportunity of suppressing the percolation of vortices by simply adjusting the coupling 
constants; we can then study the F-P observables in screened phases with and without 
percolation. We note that the "confinement-like" phase is a screened phase, rather than a 
true confinement phase, in that the energy of a static color charge is not infinite, because it 
can be screened by the dynamical Higgs particles. Nevertheless, in this phase the Coulomb 
potential is confining, and the Coulomb energy of an isolated charged state of the form 
( |2.1| ) (an unscreened charge) is infrared infinite ||. 3 

In Figs. [l(] and 11 we display the results for p(X) and F(X) on a 12 4 lattice at ft = 2.1 



and 7 = 0.6, which is inside the symmetric (or "confinement-like") phase. Data extracted 
from the full lattice, vortex-only, and vortex-removed configurations are shown on each 



3 The Coulomb potential is only an upper bound on the static quark potential [H]; a confining Coulomb 
potential is a necessary but not sufficient condition for confinement. 
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Figure 10: Eigenvalue densities for a gauge-Higgs system in the "confined", unbroken 
remnant symmetry phase. Data for the full, vortex-only, and vortex-removed configurations 
are shown, taken on a 12 4 lattice. 




Figure 11: F(X) for a gauge-Higgs system in the "confined" phase, on a 12 4 lattice. Data 
is shown for the full, vortex-only, and vortex-removed configurations. 

graph, and these are qualitatively very similar to what we have already seen in the pure- 
gauge theory (equivalent to 7 = 0). But things change drastically as we move across the 
transition into the Higgs phase. The results for p(X) and F(X) at f3 = 2.1 and 7 = 1.2, 
which is inside the Higgs phase, are shown in Figs. |l2| and Data, taken on a 12 4 lattice, 
is displayed for the full configurations only. Note that these figures look almost identical 



to the corresponding vortex-removed data in the symmetric phase, Figs. |10| and |11| which 
in turn are very close to the vortex-removed data obtained in the pure gauge theory at 
(3 = 2.1, on a 12 4 lattice. 

This last feature is worth stressing. From the point of view of the F-P operator, a 
thermalized configuration in a pure-gauge theory factors into a piece which actually does 
the confining (the vortex-only configuration) , and a piece which closely resembles the lattice 
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Figure 12: The F-P eigenvalue density in the Higgs (broken remnant symmetry) phase, for 
the full configurations. 
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Figure 13: F(\) in the Higgs phase, for full configurations. 

of a gauge theory in the Higgs phase (the vortex-removed configuration). 

Vortex-removal in the Higgs phase does not affect the data for p(A) appreciably. 4 It 
seems safe to conclude that this behavior of the F-P observables, found in the Higgs phase, 
is consistent with an infrared finite Coulomb self-energy. 

4. Eigenvalue density in the deconfined phase 

It was reported in ref. || that the instantaneous color Coulomb potential V cou i(r) is linearly 
rising at large separation, i.e. V cou i{r) ~ a cou ir where a cou i is a Coulomb string tension, 
in the high-temperature deconfined phase of pure gauge theory. This fact is surprising 
at first sight, because the free energy V(r) of a static quark pair in the deconfined phase 

4 The eigenvalue density for vortex-only configurations shows a set of sharp, narrow peaks, very much 
like the sum of delta-functions in the zero-field limit. 
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is not confining. But it is not paradoxical, because V cou i(r) is the energy of unscreened 
charges, and provides only an upper bound on V(r) Q. Confining behavior in V cou i(R) is 
a necessary but not sufficient condition for confinement, as we have already noted. 

In fact, the confining behavior of V cou i(R) in the deconfinement phase was to be ex- 
pected. The color Coulomb potential is derived from the non-local term H cou i in the 
Coulomb gauge Hamiltonian, this term depends on the Faddeev-Popov operator via the 
expression M (— V 2 )Af , and the F-P operator, on the lattice, depends only on spacelike 
links at a fixed time. But we know that even at very high temperatures, spacelike links at 
any fixed time form a confining D = 3 dimensional lattice, and spacelike Wilson loops have 
an area-law falloff just as in the zero-temperature case Jl6|]. Since V cou i(R) depends only on 
the three-dimensional lattice, it is natural that V cou i(R) confines at any temperature, and 
that remnant gauge symmetry (as opposed to center symmetry) is realized in the unbroken 
phase. 

The role of center vortices in producing an area law for spacelike loops at high tem- 
peratures has been discussed in refs. (I?]]. If the confinement property of spacelike links 
is eliminated by vortex removal, then by the previous reasoning we would expect a cou i to 
vanish as a consequence. This consequence was also verified in ref. ||. 

Given these results, it may be expected that the color-Coulomb self-energy £ r of an 
isolated static charge is infrared divergent in the deconfined phase, and that this is associ- 
ated with an enhancement of the eigenvalue density p(X) of the Faddeev-Popov operator, 
and of F(X), the expectation value of (—V 2 ) in F-P eigenstates. 

To test this we have evaluated p(A) and F(X) at j3 = 2.3 on 16 3 x 2 and 20 3 x 2 lattices, 
which is inside the deconfined phase. We have done this for the full configurations, and also 
for the vortex-only and vortex-removed configurations, defined as in the zero-temperature 
case. The results are shown in Figs. |l4] and |l^. The striking feature of these figures is their 
strong resemblance to the corresponding figures in the confined phase at zero temperature, 
namely Figs. p] and || for full configurations, Figs. [| and [| for vortex-only configurations, 
and Figs. [?] and |8] for the vortex-removed configurations. Although we have not attempted 
to determine the critical exponents in the deconfined phase, it is clear that there is an 
enhanced density of low-lying eigenvalues of the Faddeev-Popov operator, and that this 
is associated with the center- vortex content of the configuration. It is hard to avoid the 
conclusion that the Gribov scenario and the vortex-dominance theory apply in both the 
confined and deconfined phases. 

5. Thin vortices and the eigenvalue density 

In the preceding sections we have studied the eigenvalue distribution in thin vortex con- 
figurations, extracted via center projection from thermalized lattices. It is also of interest 
to study the eigenvalue spectrum of thin vortex configurations of some definite geometry, 
such as an array. 

A single thin vortex, occupying two parallel planes in the four-dimensional lattice, can 
be constructed in the following way: Begin with the zero-field configuration, i.e. all links 
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Figure 14: The F-P eigenvalue density in the deconfined phase for full, vortex-only, and 
vortex-removed configurations. 



equal to the identity matrix. Then set U2{x,y,z,t) = —I at sites 

1<^<^, y = l, aHz,t. (5.1) 

This creates two P-plaquettes in every xy plane, which extends to two vortex lines (stacks 
of P-plaquettes along a line on the dual lattice) at a fixed time, and two vortex planes 
when extended also in the i-direction of the lattice. The two planes bound a connected 
(Dirac) 3-volume, and will therefore be referred to as a single closed vortex. 

Generalizing slightly, we create N vortices parallel to the zt plane by setting 

U 2 (x,y n ,z,t) = -I (5.2) 




Figure 15: F(X) in the deconfined phase for full, vortex-only, and vortex-removed configu- 
rations. 
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Figure 16: Distribution of low-lying eigenvalues in zero-field (0,0), one-vortex (1,1), and 
vortex array (1,3), (3,3) configurations. 



for 



L 

!<*<-, 



1 + ^L, n = 0,l,...,iV-l}. 



(5.3) 



In the same way (and making a few arbitrary choices of location), we create N vortices 
parallel to the xt plane by setting 



U 3 (x,y,z n ,t) 



for 



3<y < - + 2, [zn = l + j^L, n = 0,l,...,iV-l} 
and parallel to yt plane by 



Ui(x n ,y,z,t) 



for 



< z < L - 1 



n 



1 + -L, n = 0,l,...,iV-l 



}■ 



(5.4) 

(5.5) 
(5.6) 
(5.7) 



We will consider a class of configurations characterized by the pair of integers (N, P) , where 
P is the number of planar orientations (i.e. P = 1 means only vortices parallel to the zt 
plane, while P = 3 means vortices parallel to the xt, yt and zt planes), and N is the 
number of closed vortices in each planar orientation. 

As already noted, any configuration consisting of links equal to ±7 trivially fulfills the 
Coulomb gauge condition, but also generally lies outside the Gribov region. We therefore 
perform a random gauge transformation on each thin vortex configuration, and fix to a 
Coulomb gauge copy in the Gribov region by over-relaxation. Then the F-P eigenvalue 
spectrum (the same at any time-slice) is extracted by the Arnoldi algorithm. 

The results are shown in Fig. |0] for the first 20 F-P eigenvalues obtained on a 12 4 lattice 
for zero-field (0,0), one vortex (1,1), three vortex (1,3), and nine vortex (3,3) configurations. 
As mentioned previously, it can be shown analytically that arrays of thin vortices can 
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have some additional zero modes, beyond the three trivial zero modes which exist for 
any lattice configuration. These additional modes are seen numerically in our calculation 
for the vortex spectrum. Apart from these extra zero modes, we see that the low-lying 
eigenvalue spectrum of the one- vortex configuration is not much different from the zero-field 
result. But the magnitudes of the low- lying eigenvalues changes drop abruptly, compared 
to the laplacian (0,0) eigenvalues, upon going to a three vortex configuration, and drop 
still further in a nine vortex array. This is only a qualitative result, but it does illustrate 
quite clearly the connection between vortices and the Gribov horizon. When the vortex 
geometry is chosen to imitate various features of percolating vortices, e.g. piercing planes 
in all directions and distributed throughout the lattice, then the low-lying eigenvalues have 
very small magnitudes as compared to the zero-field (or, in perturbation theory, zero-th 
order) result. This implies a pileup of F-P eigenvalues near A = 0, which we know is 
required for confinement. 

6. Convexity of FMR and Gribov regions in SU(2) lattice gauge theory 

In the section that follows this one, we will show that vortex configurations play a special 
role in the geometry of the fundamental modular region A and of the Gribov region f2. 
But first we interrupt our narrative to establish a very general convexity property of these 
regions in lattice gauge theory. 

We start by recalling the well known convexity of A and f2 in continuum gauge theory 
§. The Gribov region $1 is the set of configurations that are transverse V • A = 0, and 
for which the Faddeev-Popov operator M{A) = — V • T){A) is positive. It is also the set of 
all relative minima on all gauge orbits of the functional F^ig) = \ \ 9 A\\ 2 , where ||^4|| is the 
Hilbert norm of the configuration A, and 9 A is its gauge transform by g. The fundamental 
modular region A is the subset of the Gribov region, Acfl, which consists of all absolute 
minima Fa(1) < || 9 ^4|| 2 for all local gauge transformations g(x). It is well known that 
in continuum gauge theory both of these regions are convex namely, if A\ and A2 are 
configurations in A (or SI), then so is A = aA\ + PA2, where < a < 1, and = 1 — a. 
Geometrically, A lies on the line segment that joins A\ and A2. 

We shall establish a similar, but slightly weaker property that holds in SU(2) lattice 
gauge theory. This is quite surprising because convexity is a linear concept, whereas lattice 
configurations are elements of a non-linear space. We parametrize SU(2) configurations by 
Ui(x) = bi(x) + ia ■ Si(x), where h(x) = ±[1 — aj(x)] 1 ^ 2 . The aj(x) would be a complete set 
of coordinates, but for the sign ambiguity of bi(x) above, corresponding to northern and 
southern hemisphere. We call U+ the set of configurations U where all link variables lie on 
the northern hemisphere 

U + ^{U:b i (x) = +[l-a 2 (x)} 1 / 2 }, (6.1) 
for all x and i. The di{x) are a complete set of coordinates in U + . 
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In minimal Coulomb gauge the fundamental modular region A has the defining prop- 
erty that a configuration U in A satisfies 

J2 Re Tr [g(x)~ l Ui(x)g(x + i)] < J^Re Tr U^x) (6.2) 

XI XI 

for all g(x). Thus the gauge choice makes all the Ui(x) as close to the identity as possible, 
in an equitable way over the whole lattice. In this gauge, the link variables U{(x) for 
equilibrated configurations lie overwhelmingly in the northern hemisphere, especially in 
the continuum limit f3 — > oo. We define A + to be the restriction of A to W+, the set of 
configurations whose link variables all lie in the northern hemisphere 

A + = An^+, (6.3) 

and we call it the restricted fundamental modular region. We expect that these are the 
important configurations in the continuum limit. This is in fact necessary if the gauge-fixed 
lattice theory possesses a continuum limit. In this case 3i(x) — ► —2aAi(x), where a is the 
lattice spacing, and A{(x) is the continuum gauge connection. 

In the minimal lattice Coulomb gauge, the 3i(x) satisfy the lattice transversality con- 
dition 

[div a] b (x) ee Yl [*) " * " *)] = °- ( 6 - 4 ) 

i 

This is a linear condition on the coordinates a\(x). This suggests that we identify con- 
figurations in U + with the space of coordinates a\(x), and that we add configurations by 
adding coordinates 

a,i{x) = [aai + /?a 2 ]i(x) = aai^(x) + /3a 2 ,i(x). (6.5) 

This yields a well-defined configuration U £ IA+ only if 3%(x) < 1 for all links. This is 
assured for the case that is of interest to us where < a < 1 and (3 = 1 — a. Indeed, by 
the triangle inequality we have 

\Si(x)\ < \aai ti (x)\ + \(3a 2 ,i(x)\ <a + (3 = 1. (6.6) 

Statement: In SU(2) lattice gauge theory the restricted fundamental modular region 
is convex. More precisely if a± and a 2 lie in A+, then a = aa\ + (3a 2 also lies in A + for 
< a < 1 and (3 = 1 — a. 

The proof is given in Appendix ^. 

We may establish the same convexity property for the Gribov region f2, namely that 
O4. = flfl is convex. The Gribov region has the defining property 

^ \bi(x) [ui(x + i) - uji(x)} 2 (6.7) 

XI 

- [uji(x + i) — uji(x)] ■ di(x) x [c3i(x + t) + Ui(x)]} > 0, 

for all lo. The proof for f2 is the same as for A because this inequality has the same structure 
as ( |6.2[ ), being linear in aj(x) and bi(x). 
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7. Vortices as vertices 

We have seen that $1+ is convex, so one might think it has a simple oval shape. We 
shall show however that when thin center vortex configurations are gauge transformed 
into minimal Coulomb gauge they are mapped into points Uq on the boundary d£l of 
the Gribov region f2 where this boundary has a wedge-conical singularity of a type that 
is described below. This illustrates the intimate relation of dominance by configurations 
in the neighborhood of the Gribov horizon or in the neighborhood of thin center vortex 
configurations. 

7.1 Previous results 

We call a "thin center vortex configuration" (or center configuration) one where, in maximal 
center gauge, all link variables are center elements Ui(x) = Zi(x). [Such a configuration 
may be characterized in a gauge-invariant way by the statement that all holonomies are 
center elements. (A holonomy is a closed Wilson loop.) This definition also holds in 
continuum gauge theory] When any configuration is gauge transformed by a minimization 
procedure, such as used here, it is mapped into the Gribov region Q. In ]|] it was shown: 
When a thin center vortex configuration is gauge transformed into minimal Coulomb gauge 
it is mapped onto a configuration Uq that lies on the boundary d£l of 0. Moreover its 
Faddeev-Popov operator M(Uq) = —V -T>{Uq) has a non-trivial null space that is (N 2 — 1)- 
dimensional. (This is in addition to the (iV 2 — l)-dimensional trivial null-space consisting 
of constant eigenvectors. Here and below we do not count trivial eigenvectors that are 
generators of global gauge transformations.) Likewise, when an abelian configuration is 
gauge-transformed into the Gribov region, its Faddeev-Popov operator has a non-trivial 
i?-dimensional null space, where R is the rank of the group. The reason is that the gauge 
orbit of a center or abelian configuration is degenerate so the gauge-invariant equation 
T>i(Uo)uj = 0, which holds for i = 1,2,3 simultaneously, has N 2 — 1 or R non-trivial 
solutions respectively. 

7.2 Tangent plane to Gribov horizon at a regular point 

Let a\{x) be a set of coordinates of the group element Ui(x) = U[a\{x)]. It is convenient 
to use coordinates a\(x) such that the Coulomb gauge condition is linear, (V • a) b (x) = 
J2i l a i( x ) ~ a i( x ~ *)] = 0, and in terms of which Q + (see above) is convex. We write 
Uo+SU = U (ao+Sa). Here 5a is an arbitrary (transverse) small variation of the coordinates 
at oq. Such a variation is a tangent vector at ao, and the space of tangent vectors constitutes 
the tangent space at ao- 

Let Uq be a configuration in Coulomb gauge, that lies on the boundary dCl of the 
Gribov region. By definition, the corresponding Faddeev-Popov operator has a non-trivial 
null eigenvector 



all other eigenvalues being non-negative. We are interested in the case of a center config- 
uration in Coulomb gauge where the non-trivial null eigenvalue of M(Uq) is (iV 2 — l)-fold 
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degenerate, but for orientation we first consider the case where the non-trivial null eigen- 
value is non-degenerate. Let Uq + 8U be a neighboring point that is also on the Gribov 
horizon, so M(Uq + SU) also possesses a null vector 

M{U + 5U)(lu + 5uj) = 0. (7.2) 

We wish to find the condition on 5a that holds whenever U (ao + 5a) also lies on the 
Gribov horizon dQ. We have 



where 



and 



M(U + 5U) = M(U ) + 5M, 
5M = Y J 5o b l {x) d xib M , 



xib 



- dM 
o X ib M o = b 



a=ao 



To first order in small quantities we have 

5Mlo + M(U )5lo = 0. 
We contract this equation with too and obtain 

{coo, 5M loq) = 0, 



(7.3) 
(7.4) 

(7.5) 
(7.6) 
(7.7) 



or 

J2 ^(x) (LO ,d xib M lo ) = 0. (7.8) 

xib 

Geometrically, this is the statement that 5a\{x) is perpendicular to the vector (u>q, d x n,M§ loq), 
so this vector is the normal to d£l at Uq = U{oq). It defines the a hyperplane that is tangent 
to the Gribov horizon at ao- 

7.3 Center vortices as singularities of the Gribov horizon 
7.3.1 General idea 

Suppose that Uq = U(ao) is a point on the Gribov horizon and that the null eigenvalue of 
M(ao) is P-fold degenerate, 

M(a ) 4 n) = 0; n = 1, . . . , P. (7.9) 

As noted above, this happens for every gauge copy in Coulomb gauge of a thin center 
vortex configuration, where P = iV 2 — 1. Under the small perturbation M(ao + 5a) = 
M(ao) + 5M , where 5M is given above, the P-fold degenerate null-space splits into P levels 
with eigenvalues <5A n , for n = 1, . . . , P, that depend linearly, 5\ n = Ylxia i( x )& a \ ( x )> on 
the 5a. 
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For a point ao on the boundary <9fi we define the Gribov region of the tangent space 
at ao to be the set of tangent vectors that point inside Q. We designate it by O ao . More 
formally, it is the set of tangent vectors 5a at ao such that ao + 5a lies in f2, 

Oa = {5a :a + 5a£ O}. (7.10) 

Recall that by definition the Gribov region Q consists of (transverse) configurations U(a) 
such that all eigenvalues of M(a) are positive. It follows that £l ao is the set of 5a such that 
all P eigenvalues 5X n (5a) are positive, 

n ao = {5a : 5X n {5a) > for n=l,...,P}. (7.11) 

This condition is quite restrictive because, for generic 5a, some number v of eigenvalues 
5X n are negative while P — v are positive, where v = 0, 1, . . . , P. As a result the boundary 
dft ao of the Gribov region at ao is not simply a tangent plane through ao, as before. Rather 
it is a high dimensional wedge-conical vertex whose shape we shall find. 

7.3.2 Degenerate perturbation theory 

The eigenvalues in an infinitesimal neighborhood of a point Uq on the Gribov horizon are 
determined by the eigenvalue equation 

[ M(U ) + 5M] ^ Y c «4 n) + = sx ( £ c -4 n) + &j j , (7.12) 



or 



oM^c„cjJ n) + M(U )5lo =5X ^c n u;J n) - ( 7 - 13 ) 



Upon contracting this equation with lo^, we obtain, by ([f!s|), the P-dimensional eigen- 
value equation 



y^fow c n = 5X Cm, (7.14) 



where 



5a mn = ( oJq" 1 ^ i ^ w o™' 1 



xib 



5a\(x) (4 m) ,^ i6 M Q 4 n) ). (7.15) 



Abstractly, <5a is a tangent vector in lattice configuration space at point ao in lattice 
configuration space; the components of this vector, in a suitable basis, are denoted 5a\(x). 
We may also regard the set of numbers 5a mn as components of 5a in some other basis. 

The eigenvalue equation has P solutions 5X n . The condition that they all be positive, 
which determines f2 ao , the Gribov region at ao, is the condition that the matrix 5a mn 
define a strictly positive form. Such a form is characterized by the Sylvester criterion that 
det 5a mn be positive, together with the determinant of all principle minors (diagonal square 
submatrices) of this matrix. The boundary of this region is determined by the condition 
that one eigenvalue vanish, which happens when the determinant vanishes, 

det 5a mn = 0. (7-16) 
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7.3.3 Two-fold degeneracy 

We first analyze 2-fold degeneracy. In this case positivity of the two principle minors and 
of the determinant reads 

5a n > 0, 5a 2 2 > 0, Sa u 5a 2 2 - 5a\ 2 > 0. (7.17) 

In terms of 5a + = \{5a\\ + 6022) and <5a_ = \{5a\\ — 60,22), the last condition reads 

5a\ - 5a 2 _ - 5a\ 2 > 0. (7.18) 

The three inequalities define the interior of the "future" cone in the 3- variables 5a + , 5a- 
and 5a\2, with vertex at the origin. Thus the boundary of the Gribov region at ao is a 
cone in these 3 variables. Taking account of the remaining components of 5a, the conical 
singularity can be viewed as a kind of wedge in higher dimensions. 

7.3.4 Three-fold degeneracy 

For SU(2) gauge theory the thin-vortex configurations are 3-fold degenerate points on the 
Gribov horizon. In this case the positivity of the principle minors and the determinants 
reads 

5an > 0, 5a 22 > 0, 5a 33 > (7.19) 



San 5a22 — 5a\ 2 > 
5a\\ 5a%3 — 5a\^ > 

5a 2 2 5a 33 - 5al 3 > (7.20) 

detoa mn >0. (7.21) 

These 7 inequalities characterize the Gribov region Q a() of the tangent space at a point ao 
on the Gribov horizon that is 3-fold degenerate. 

From our discussion of 2-fold degeneracy we know that the positivity of each of the 
2 by 2 determinants and of its diagonal elements defines a "future" cone. We call these 
future cones F12, F13 and F23 respectively. For example F12 is defined by San > 0, and 
5a22 > 0, and 

[\{5 ail + 5a 2 2)] 2 - [k(6a u - 5a 2 2)} 2 - 5a 2 12 > 0, (7.22) 

etc. Each is a cone of opening half-angle 7r/4, and we have shown that O ao is contained in 
the intersection of the 3 future cones, 

O.,, C F 12 nF 13 nF 23 . (7.23) 

This condition and the 3 by 3 determinantal inequality det 5a mn > characterize the 
Gribov region Q ao in the tangent space of a point Uq = U (ao) that is 3-fold degenerate. 
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7.3.5 Over-all picture of the Gribov horizon and its center-vortex singularities 

We have seen that J)+ is convex and that thin center vortex configurations are wedge- 
conical singularities on the boundary of d£l. Those that are on are extremal elements, 
like the tips on a very high-dimensional pineapple. Indeed, each center configuration is an 
isolated point. If one moves a small distance from a center configuration it is no longer a 
center configuration. Its gauge transform U(ao) in Coulomb gauge is likewise an isolated 
point on the Gribov horizon. Thus the wedge in the boundary d£l at ao that we have 
just described occurs at an isolated point where the Gribov horizon may be said to have a 
"pinch". 

In SU(2) gauge theory there are 2 dV center configurations (where d is the number 
of dimensions of space, and V is the volume of the lattice) because there are dV links 
in the lattice and there are 2 center elements in SU(2). These are related by 2 V gauge 
transformations, so there are 2^ d ~ l ^ v center orbits. The absolute minimum of each of these 
orbits lies on the common boundary of the fundamental modular region dA and the Gribov 
region d£l. So there are at least 2^ d ~^ v tips on the above-mentioned pineapple, a truly 
enormous number. Moreover for each such orbit there are many Gribov copies, all lying 
on dO, (spin glass problem). These are all singular points of the Gribov horizon of the type 
described. For SU(2) there may not be any other singular points on d£l. It is possible that 
the thin center vortex configurations provide a rather fine triangulation of d£l. 

Note: The gauge transform in Coulomb gauge of an abelian configuration also lies on 
the Gribov horizon. However the SU(2) group is of rank 1, so, for SU(2) gauge theory, 
such a configuration is invariant under a one-parameter group of transformations only. As a 
result, the corresponding null space M(Uq)ujq = is only one-dimensional, and the present 
considerations do not indicate that these are singularities of d£l. 



8. Coulomb gauge as an attractive fixed point of the interpolating gauge 

The data reported above have been obtained by numerically gauge fixing to the minimal 
Coulomb gauge using a well-defined numerical minimization procedure on a lattice of finite 
volume. However not every lattice gauge corresponds to a well-defined continuum gauge, 
so one may well ask whether our results have a continuum analog, especially since in 
continuum theory the Coulomb gauge is a singular gauge. Indeed, in Coulomb gauge, some 
propagators of the form 1/k 2 appear instead of l/(fcQ+k 2 ), which leads to unrenormalizable 
divergences J dko in some closed loops. This has been treated by using the action in phase 
space (first order formalism), which allows a systematic cancellation of these divergences 



1 18], or by using an interpolating gauge [19] with gauge condition 



a<9 A) + V • A = 0. (8.1) 

This gauge interpolates between the Landau gauge a = 1 and the Coulomb gauge, a = 0, 
and may also be achieved by a numerical minimization procedure. For a > it is a regular 
continuum gauge, and the gauge parameter a provides a regularization of the divergences 
of the Coulomb gauge, which is obtained in the limit a — > of the interpolating gauge. 
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However as a possible obstacle to this regularization, it must be noted that the quan- 
tities that appear in the interpolating gauge condition ( |8.1| ) are unrenormalized quantities, 
and the unrenormalized gauge parameter depends on the ultraviolet cut-off A, 

a = a(A/ Aqcd)- (8.2) 

For the Coulomb-gauge limit of the interpolating gauge to be a smooth limit it is necessary 
that the gauge parameter a (A/ Aqcd) flows toward (or at least not away from) the Coulomb- 
gauge value a = 0, as the ultraviolet cut-off is removed, A — » oo. In Appendix [C] its 
dependence on A is calculated in the neighborhood of a = to one-loop order in the 
perturbative renormalization group, with the result in pure SU(A") gauge theory, 

.... , const 

° (AMqcd) = HA/WAI ' < 8 ' 3 ' 



This gives 



lim o(A/A QC d) = 0, (8.4) 

A— >oo 



and we conclude that in some neighborhood of a = 0, the Coulomb gauge value o = 
is an attractive fixed point of the renormalization-group flow. This removes the possible 
obstacle to regularization of the Coulomb gauge by the interpolating gauge. 

[Similarly the renormalized gauge parameter a R = a R (fi/ Aqcd), depends on the renor- 
malization mass /u. As the renormalization mass [i gets large, which is the appropriate 
choice for calculations at high momentum, one has at one-loop order 

oV/Aqod> = wCr 1 (8 ' 5) 

So again, the Coulomb-gauge value a R = is an attractive fixed point of the renormalization- 
group flow.] 



9. Conclusions 



We have found that the low-lying eigenvalues of the Faddeev-Popov operator, in thermal- 
ized lattices, tend towards zero as the lattice volume increases. This means that in the 
infinite volume limit, thermalized configurations lie on the Gribov horizon. That fact alone 
would not allow us to make any strong conclusions about the energy of unscreened color 
charge. However, the data also indicate that the density p(X) of F-P eigenvalues goes as 
a small power of A, at infinite volume, as A — ► 0. Together with the behavior of -F(A) at 
A —* 0, we conclude that the energy of an unscreened color charge is infrared divergent, 
and that this divergence can be attributed to the near-zero modes of the Faddeev-Popov 
operator. 

This evidence clearly supports the Gribov horizon scenario advocated by Gribov and 
Zwanziger in ref. This scenario was invented to account for confinement, and the 

reader may be surprised to find a linearly rising color-Coulomb potential in the deconfined 
phase. However this is nicely explained by the horizon scenario in Coulomb gauge, as 
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we now explain. In Coulomb gauge the gauge fixing is done independently on each 3- 
dimensional time slice, and may be done on a single time slice. According to the horizon 
scenario, on each time slice, 3-dimensional configurations Ai(x) are favored that lie near 
the horizon of a 3-dimensional gauge theory, and this enhances the instantaneous color- 
Coulomb potential. This is true for every temperature T, including in the deconfined phase, 
because temperature determines the extent of the lattice in the fourth dimension. Thus, 
the horizon scenario provides a framework in which confinement may be understood, but it 
is not detailed enough to tell us under what conditions the infinite color-Coulomb potential 
may be screened to give a finite self-energy, as measured by the Polyakov loop. 

By factoring thermalized lattices into vortex-only and vortex-removed components, we 
have also been able to show that the constant density of low-lying F-P eigenvalues can be 
entirely attributed to the vortex component. We find that the eigenvalue density of the 
vortex component is qualitatively similar to that of the full configuration. The eigenvalue 
density of the vortex-removed component, on the other hand, is drastically different from 
that of the full configuration. This density can be interpreted as simply a small perturbation 
of the zero-field (or zero-th order) result, and it is identical in form to the (non-confining) 
eigenvalue density of lattice configurations in the Higgs phase of a gauge-Higgs theory. 

These findings establish a firm connection between the center vortex and the Gribov 
horizon confinement scenarios. According to the center vortex doctrine, fluctuations in 
the vortex linking number are responsible for the area law falloff of Wilson loops. It now 
appears that vortex configurations are also responsible for the enhanced density of near-zero 
F-P eigenvalues, which is essential to the Gribov horizon picture. This result is consistent 
with previous results in ref. [20], where it was found that vortex removal also removes the 
Coulomb string tension of the color Coulomb potential. It is also consistent with recent 
investigations of Gattnar, Langfeld, and Reinhardt in Landau gauge. 

The F-P eigenvalue spectrum at high temperatures, with and without vortices, turns 
out quite similar to the corresponding results at low-temperature. This similarity was to be 
expected, since the F-P operator depends only on spacelike links at a fixed time, and even 
at high T these form a confining three-dimensional ensemble for spacelike Wilson loops. 
The Gribov scenario must therefore be relevant to physics in the deconfined phase; cf. ref. 
1 22] for a recent application. 

We also report a result which supports the consistency of Coulomb gauge itself in the 
continuum limit. Coulomb gauge is very singular in continuum perturbation theory, and 
one method of making it better defined is to view Coulomb gauge as a non-singular limit of 
the more general gauge condition ad^A^ + V • A = 0. The success of this approach depends 
on whether the Coulomb gauge limit, a = 0, is an attractive ultraviolet fixed point of the 
renormalization group flow. Here we have shown that this requirement is satisfied. 

Finally, we have uncovered an intriguing geometrical property of thin vortices in lattice 
configuration space. In ref. || it was shown that thin vortices (gauge equivalent to center 
configurations) lie on the Gribov horizon. The Gribov horizon is a convex manifold, and we 
have shown here that thin vortices are conical singularities on that manifold. Percolating 
thick vortices appear to be ubiquitous in thermalized lattice configurations at or near the 
Gribov horizon; it is conceivable that the special geometrical status of thin vortices is in 
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some way related to the ubiquity of their finite-thickness cousins. 
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A. Finite-volume scaling of low-lying F-P eigenvalues 

In N x N random matrix models, the tail of the eigenvalue distribution often displays a 
universal scaling behavior with N. This fact has found important applications in the study 
of chiral symmetry breaking, where this sort of universal scaling is found in the density of 
near-zero eigenvalues of the Euclidean Dirac operator as a function of lattice volume (which 



is proportional to the number of eigenvalues) [23]. In our case, we are also interested in the 



density of near-zero eigenvalues in the infinite 3- volume limit. The statement in this case 
is as follows (cf., e.g., ref. f24|): Suppose the normalized density of low-lying eigenvalues, 
at very large volumes, goes as 

p(A) = K \ a (A.l) 

where k is some constant. Then the density of eigenvalues, the average spacing between the 
low-lying eigenvalues, and the probability density P(A n ) of the n-th low-lying eigenvalue, 
agree for every lattice 3- volume V3 , if the eigenvalues themselves are rescaled according to 

z = \V.^. (A.2) 

The argument goes as follows: The number of eigenvalues N[X, AA] in the interval [A — 
±AA, A + |AA] is 

N[X,AX] = 3V 3 p(X)AX. (A.3) 
Then in terms of rescaled eigenvalues z = XV£, we have 

N[X, AA] = 3KV^ p(1+a) z a Az. (AA) 

If we require that this number of eigenvalues depends only on the rescaled variables z, Az, 
then it is necessary that p = 1/(1 + a), and eq. ( |A.2j ) follows. 



Our strategy is to plot the probability density P(\ m i n ), rescaled by a factor V 3 l / i - l+a \ 
as a function of the variable z m i n = X m i n V^'^ +a \ where \ m in is the lowest non-zero 
eigenvalue. The rescaling of the probability density ensures that its integral over z m i n is 
unity. P(Xmin) is computed on 8 4 to 20 4 lattice volumes, at various values of a. If we 
can find a value of a for which the (rescaled) probability densities P(X m i n ) coincide at all 
lattice sizes, then this determines a, and in turn the behavior of the eigenvalue density 
near A = 0. 
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The results for rescaled P(X m in), for both the full configurations {\ m in = A4), and 
the vortex-only configurations (A m j„ = A7), are displayed in Fig. 17. We show only three 
values of a in each case, which include our best estimate for the scaling a. We find 
that a = 0.25 ± 0.05 for the full configurations, and a = 0.0 ± 0.05 for the vortex-only 
configurations seem to give the best scaling results; the error estimate is subjective. 

The lattice data for L = 20 in the vortex-only configurations calls for some comment. 
This distribution does not match the distributions obtained at smaller L at any a. This 
could mean that the scaling hypothesis for low-lying eigenvalues is invalid, but in our 
opinion the mismatch at L = 20 has a different explanation: We believe the problem is 
connected to difficulties we have encountered with gauge fixing on these large lattices. As 
lattice size is increased, typical values of \ m in become smaller, and the number of over- 
relaxation steps required to satisfy our Coulomb gauge convergence criterion increases. In 
general, on a given lattice size, the number of over-relaxation steps required to gauge fix 
is inversely correlated with the size of \ m in in the gauge-fixed configuration. 

It happens occasionally that even after 10,000 gauge-fixing steps on a given time-slice, 
the spatial lattice configuration has not converged to Coulomb gauge. When this happens, 
we perform a random gauge transformation at that time slice and gauge fix a second time, 
to a different Gribov copy (note that in Coulomb gauge, time-slices can be, and were, gauge- 
fixed independently). This procedure most likely biases the result towards higher average 
values of \ m in- In doing things this way, we are almost certainly modifying the probability 
distribution in Gribov region, giving lower measure to Gribov copies which are closer to 
the horizon. On the other hand, simply excluding these hard-to-gauge-fix lattices would 
probably introduce an even worse bias. On smaller (L < 20) vortex-only lattices, and on 
the full, unmodified lattices, problems with convergence to Coulomb gauge are uncommon. 
However, on the L = 20 vortex-only lattice, there is a convergence failure on almost 38% of 
all time slices; on these slices we have performed a random gauge transformation to move 
to a different Gribov copy. The rate of convergence failure on vortex-only 16 4 lattices is 
five times lower, and the rate of failure on unmodified 20 4 lattices is eight times lower, 
than on the vortex-only 20 4 lattice. For this reason, we believe that the bias towards larger 
eigenvalues is by far the worst on the vortex-only 20 4 lattices, and this is in fact where we 
see the mismatch, in P(A m j n ) vs. z m i n at a ~ 0, with the other vortex-only lattice volumes. 

Note that for the vortex-only configurations, a = is consistent with a finite, non-zero 
value for p(0), as shown in eq. ( p. 25 ). On the other hand, it is not excluded that a could in 
fact be slightly negative, in which case p(X) actually diverges at A = 0. Since the eigenvalue 
density in Fig. || does appear to actually rise as A — > 0, before suddenly falling, a divergent 
behavior at A = in the infinite volume limit is not at all excluded. 

Next we consider F(X) as A — » 0. We have fit the average value of F at the lowest 
non-zero eigenvalue, (F(X m i n )), to the form 

(F(X min )) = ^ + b. (A.5) 

Our result, within errorbars, is consistent with 

(-F(A4)) = — full configurations, 

1j 
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Rescaled P(A, min ) 

Full configurations Vortex-only configurations 




Figure 17: Probability distribution of the values of the lowest non-trivial FP eigenvalue 
^min on a variety of lattice sizes 8 4 — 20 4 , as a function of the rescaled eigenvalue z m i n = 
\ m i n L 3 / 1+a . Results for the full configurations are in the left-hand column of figures, with 
the vortex-only results in the right-hand column. Each figure corresponds to a different 
choice of a; a = 0.1,0.25.0.4 for the full configurations, a = —0.2,0.0,0.2 for the vortex- 
only configurations. The closest match of rescaled P{\ m in) at different volumes comes at 
a k, 0.25 for the full configurations, and a k, for vortex only configurations (with the 
exception of L = 20, see text). 
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Figure 18: F(X) as in Fig. |!| together with a best fit to eq. ( |A.8| ). 



20 vortex-only 
(F(Xj)) = 1 + — j-j configurations. 

Li ' 



(A.6) 



Since A m j„ — ► at infinite volume, then jP(0) = 1 for vortex-only configurations, as stated 
in eq. ( |2.25| ). For the full configurations, since we estimate a ~ 0.25, with the consequence 
that (X m i n ) ~ 1/L 2 - 4 , it is reasonable to guess that for the full configurations near A = 0, 



F(X) ~ A ' 42 . 
We have therefore tried the following fit 

F(X) = aX p + bX 



(A.7) 
(A.8) 



to the L = 20 data, where the linear term is motivated from the fact that perturbative 
behavior (linear dependence on A ) is expected at large A . The fit gives the result shown 
in Fig. [ll| with an exponent p = 0.38 which is not far off our guess of p = 0.42. 

B. Proof of convexity of FMR in SU(2) lattice gauge theory 

The defining property ( |6.2[ ) of the fundamental modular region may be written 

Re Tr [Ui(x) gi (x)] < £ Re Tr U^x), (B.l) 

xi xi 

where g%{x) = g(x + i)g(x) . We write gi(x) = Ci(x) + ia ■ g*i(x), where Ci{x) = ±[1 — 
•^(x) 2 ] 1 / 2 , and the defining property reads, 

-^2ai(x)-gi(x) <J2Hx)[t-d(x)}, ( B - 2 ) 

xi ix 

where the last factor is positive, 1 — q(x) > 0. Given that this property holds for configu- 
rations a\ and 02, 



ai,i(x) ■ gi(x) < bi t i(x) [1 - Ci(x)] 



(B.3) 
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- &A X ) ■ 9t(x) < ^ b 2ii (x) [1 - a{x)} , (B.4) 

xi xi 

then we will show that it holds for a = aa\ + (3a 2 , as stated in section V. Upon multiplying 
the first equation by a and the second by (3 and adding we obtain, 

- a,i(x) ■ gi(x) < [ab lt i(x) + (3b 2)i (x)\ [1 - , (B.5) 



where 3{(x) = aai i j(x) + /3a2 i j(x). It follows from this that the defining property ( |B.2| ) that 
we wish to establish will be proven if we can show that on each link (xi) the inequality 

q6 m (x) +Pb 2 ,i{x) < bi{x) (B.6) 

holds, where h(x) = [1 - a^x) 2 ] 1 / 2 . 

We introduce trigonometric functions 

h,i( x ) = cos 6i,ii x )'' ai,i( x ) = sm8^i(x) h lti (x), (B.7) 
h,i( x ) = cos ^2,i( x )i a 2 ,i(x) = sin6>2,i(x) n 2>i (x), (B.8) 
= cos 9i(x); a,i(x) = sm9i{x) fii(x), (B.9) 

where < 9 lti (x) < vr/2, and < 9 2 ,i(x) < it/2, and < 9i(x) < vr/2. They are related by 

asin#j(x) n\^(x) + (3 sin 9 2 ^(x) n 2 ^(x) = sm9i(x)hi(x), (B.10) 

and the property to be proved reads 

acos9\^(x) + 13 cos 9 2 ,i (x) < cos#j(x). (B.ll) 

We have 

q cos 9i t i(x) + cos 9 2t i(x) < cos [a9i t i(x) + I39 2 ^(x)] , (B.12) 
by the convexity of the cosine function, so it is sufficient to prove 

cos [aBx,i{x) + P&2,i(x)\ < cos9i(x). (B.13) 



The triangle inequality applied to (|B,10D reads 



s'm9i(x) < a sin 9\^(x) + /3sin 9 2 ^(x) 

< sm[a9 1 4x)+(39 2 , l (x)}, (B.14) 

where we have also used the convexity of the sine function. This yields the inequality 
9%{x) < a9± : i(x) + (39 2: i(x), from which ( |B.11| ) follows. This proves the convexity of A + . 
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C. Renormalization-group flows toward Coulomb gauge 

We compute the flow 5 of the gauge parameter a that appears in the interpolating gauge 
condition ( |8.1| ) in the neighborhood of a = 0. 

Lorentz-invariance is not manifest in the interpolating gauge, and Aq and A renormal- 
ize independently, 

A = Z Ao A R ; A = Z A A R , (C.l) 

where superscript R designates renormalized quantities. The renormalized gauge parameter 
a R is defined by 

a = Z a a R , (C.2) 
where the renormalization constant Z a is determined by the renormalized gauge condition, 

a R d A R + V • A R = 0. 

This gives 

Z a = Z^Z A . (C.3) 

The 7-functions are defined by 

d t lnZ a = 7 a (5), 

d t lnZ Ao = lAo (g), (C.4) 
d t lnZ A = 7a(s0, 

where t = In A, and A is the ultraviolet cut-off. The derivative is taken at fixed renormalized 
coupling constant g R . We have 

8 t In Z a = -8 t In Z Ao + 8 t In Z A , (C.5) 

which gives 

Ja(g) = -JA (g) + 7A(ff)- (C.6) 
From the perturbative expansions 

Z a = l + c a lnA (g R ) 2 + ..., 

Z Ao = l + c Ao lnA(g R ) 2 + ..., (C.7) 
Z A = 1 + caIuA (g R ) 2 + ..., 

we obtain 

la = c a g 2 + . . . , 

7A = c Ao g 2 + ..., (C.8) 
7A = c A g 2 + ■ • • , 



The calculation reported here was motivated by correspondence with Kurt Langfeld and Laurent Moy- 
aerts. 
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and 



Ca — ~ c Aq + ca- 
To one-loop order the renormalization constant satisfies 

d t In Z a = c a g + • . 



2b t 



+ ..., 



(C.9) 



(CIO) 



where &o is t ne first coefficient of the /3-function. In pure SU(iV) gauge theory it has the 
value 

*» = ^. (en) 



3 16vr 2 ' 



We obtain, to one-loop order, 



Z a = const t 



c a /2b 



(C.12) 



which, with t = In A, gives the leading dependence of the bare gauge parameter on the 
cut-off A, 

a(A/A QC D) = const (A/A Q c D ) Ca/26 °. (C.13) 

Clearly a = is a fixed point of the renormalization group. To see if it is also a 
stable fixed point, we may evaluate the coefficient c a at a = 0, namely in Coulomb gauge, 
assuming that c a is smooth in the neighborhood of a = 0. The renormalization constants 
Za and Zj± in Coulomb gauge, are given to one-loop order in eq. (B.37) of [18]. The 
coefficients have the values 



CA 



11 N 



1 N 
CA ~ 28^- 



which gives 



and 
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